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We investigate the relation between the cooperative length and the relaxation time, represented re¬ 
spectively by the culling time and the persistence time, in the Fredrickson-Andersen, Kob-Andersen 
and spiral kinetically-constrained models. By mapping the dynamics to diffusion of defects, we find 
a relation between the persistence time, Tp, which is the time until a particle moves for the first 
time, and the culling time, Tc, which is the minimal number of particles that need to move before 
a specific particle can move, Tp = r^, where 7 is model- and dimension dependent. We also show 
that the persistence function in the Kob-Andersen and Fredrickson-Andersen models decays subex- 
ponentially in time, P{t) = exp |^— (t/T)^j , but unlike previous works we find that the exponent ^ 
appears to decay to 0 as the particle density approaches 1. 
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I. INTRODUCTION 


Increasing the density of particles in granular materi¬ 
als causes them to undergo a transition from a fluid-like 
state, in which the particles can move relatively freely, 
to a jammed state, in which almost none of the particles 
can move Ea. In glasses, a similar transition occurs 
when the temperature is decreased [sl, 0] • As the mate¬ 
rial nears the glass or jamming transition, the system’s 
relaxation time increases dramatically, until it diverges 
at the critical point [H . 

The various kinetically-constrained models lili3 cap¬ 
ture the essence of the glass or jamming transitions and 
there has been much recent activity on them. Some of 
these models simulate the way that particles block each 
other’s movement by requiring that a particle can move 
only if its neighbors satisfy some condition [nMii. Other 
models add driving forces which simulate the resistance 
of jammed systems to external forces [iO-il- In gen¬ 
eral, the system is coarse-grained to a lattice, and each 
site is in one of two states, 0 or 1. In lattice-gas models 
a site in state 1 represents a particle which may move 
to an adjacent vacant site, represented by state 0, if its 
local neighborhood satisfies some model-dependent rule. 
In spin-facilitated models state 1 represents a high den¬ 
sity region in granular systems and an inactive region 
in glasses, while state 0 represents either a low density 
region or an active region in granular matter and glass¬ 
forming liquids respectively. A site can change its state 
from 0 to 1 and vice versa, with a temperature-dependent 
rate if the site’s local neighborhood satisfies some model- 
dependent rule. 

In this paper we consider the Kob-Andersen (KA) 
[m and Fredrickson-Andersen (FA) 11,123 kinetically- 
constrained models on one- and two-dimensional square 
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lattices. In the FA spin-facilitated model, a site can 
change its state from 0 to 1 and vice versa if it has at least 
m neighboring vacancies. In the KA lattice-gas model, a 
particle needs at least m adjacent vacancies before and 
after the move in order to move to a nearest neighbor 
vacant site. Higher dimensional versions of these mod¬ 
els with higher values of m have also been investigated 


The glass or jamming transitions result from coopera¬ 
tive dynamics, in the manner that particles are blocked 
by their neighbors, which in turn are blocked by their 
neighbors, and so on, such that in order for a single par¬ 
ticle to move, many others need to move before it. The 
number of these “shells” and their weight, represent the 
structural changes in the system as it nears the critical 
point, and they diverge at the critical point. In effect, 
they represent the minimal number of steps needed for 
a particle to move, which may be found by culling the 
shells iteratively. Above the critical density, or equiva¬ 
lently below the critical temperature, some of these shells 
cannot be culled since the particles in them block each 
other. This culling process is the usual manner to check 
whether a system is jammed or not, because if no shells 
remain after the culling then all the particles may move 
and the system is not jammed. The culling time repre¬ 
sents a length scale related to relaxation of the system 
dill- We note here that this length scale is not the 
only way to quantify the relation between the structure 
of the system and its dynamics jsol-l^. and it remains 
an open question which structural order parameter is a 
better choice. 

In most previous works regarding the FA and KA mod¬ 
els, the relaxation time was measured by the two-time 
density autocorrelation function {sl-lioj. In this paper 
we use the persistence function, defined as the fraction 
of particles that have not yet moved until time t (in 
lattice-gas models) or the fraction of sites that have not 
changed state until time t (in spin-facilitated models). 
The persistence function was thoroughly investigated in 
the relatively simple m = 1 models |4ll445j. but there 



are also works on higher values of m [Mil , and other 
kinetically-constrained models [22|. Generally, the den¬ 
sity autocorrelation function and the persistence function 
behave similarly. 

In this paper we study the relation between the culling 
time and the relaxation time, obtained from the persis¬ 
tence function, and show that near the critical point the 
relation is a model-dependent power law which can be ex¬ 
plained as a diffusion of rare droplets. We show that this 
is a general result by also considering another kinetically- 
constrained model, the spiral model In Section 

[Uwe describe the models investigated in this paper. Our 
results for the culling time and the relaxation time are 
shown in Sections uni and 11^ respectively, and are com¬ 
pared in Section O Section IVll summarizes the paper. 


II. THE MODELS 


We consider the KA and FA models on a d-dimensional 
square lattice. At time t = 0, each site in the lattice is 
either in state 1 with probability p, or in state 0 with 
probability v = 1 — p without correlations between sites. 
In this way, we probe the equilibrium distribution of the 
system. In the FA model, a site can change its state from 
0 to 1 and vice versa if it has at least m neighboring 
vacancies. In the KA model, sites at state 1 are occupied 
by particles, and sites at state 0 are vacant. A particle 
needs at least m adjacent vacancies before and after the 
move in order to move to a nearest neighbor vacant site. 
We consider here three cases: d = m = 1, d = 2,m = 1, 
and d = m = 2. The first two cases (m = 1) in the KA 
model are equivalent to the simple symmetric exclusion 
principle (SSEP) model [H, in which a particle can 
move if it has a neighboring vacancy. When m = 1 all 
the particles are able to move eventually (in the KA mod 
el) or change their state eventually (in the FA model), 
while if m = 2 there is a system-size-dependent value of 
the density above which a finite fraction of the particles 
will not be able to move (KA model) or change their 
state (FA model) [26|. In square systems of size L x L, 
this critical vacancy density is given by [52| 


Vc = 


A 


( 1 ) 


where A is a weak function of L [El, [EJ • In the system 
sizes we consider here, A ~ 0.25, whereas in the limit 
1/ ^ oo it is equal to 7r^/18 ~ 0.55. 

We perform on these systems two types of dynamics: 
culling dynamics and real dynamics. In the culling dy¬ 
namics, we iteratively remove the particles which are able 
to move (KA), or change to 0 the state of the sites which 
are able to do so (FA). In the real dynamics, every time 
step dt = 1/N, with N being the number of sites, one of 
the sites is chosen randomly. 

In the FA model, in order for a site to change its state, 
it first must have m neighboring vacancies as noted be¬ 
fore. If this condition is satisfied, the site changes its state 


from 0 to 1 with probability Wqi and from 1 to 0 with 
probability Wiq. In order to maintain detailed balance 
while maximizing the transition probabilities, we set 

Woi = min (l, 0 , 

Wio = min ^1, ^ . (2) 

In the KA model, if the chosen site is occupied, a random 
direction is also chosen, and the chosen particle can move 
in that direction if the neighboring site in that direction 
is empty, and the particle has at least m neighboring 
vacancies before and after the move. 

In the real dynamics, we use a continuous time, or 
rejection-free algorithm since at high densities the prob¬ 
ability that an allowable move is randomly chosen is very 
small. In this algorithm we randomly generate the num¬ 
ber of time steps that have passed between successive 
moves based on the probability that a move is possible. 
In this way, we do not wait for long periods of time until 
a move is made, but rather advance the clock in large 
random steps. 

For the culling dynamics we define the culling time cu¬ 
mulative distribution as the fraction of sites 

that started in state 1 and didn’t change to 0 until itera¬ 
tion s of the culling process, and for the real dynamics we 
define the persistence function p(^’’^)(t) as the fraction 
of particles that have not yet moved (KA) or the fraction 
of sites that started from state 1 and did not change to 
0 (FA) until time t. Obviously M(0) = P(0) = 1 for all 
models. The culling time and the persistence time 

^ respectively as the average number of it¬ 

erations needed to cull a particle and the average time 
until a particle moves (KA) or a site changes its state 
(FA) for the first time, are given by 


(d,m) 

' c 


(d,m) 

‘p 


oo 


E 


1 - M^d,m) 



_p(d,m) 
1 - Pid,m) (oc) 


( 3 ) 


where = 0 if the system is un¬ 

jammed, i.e. that all of the sites (particles) will be able to 
flip (move) eventually, and P^^’^^(oo) > > 0 

if the system is jammed, i.e. that some of the sites 
(particles) will never be able to flip (move). Therefore, 
P(^’’^)(oo) and act as the system’s Edwards- 

Anderson order parameter [E^. In the EA models, flip¬ 
ping sites only to 0 may occur in the real dynamics, al¬ 
beit with a negligible probability, and thus = 

P(^’’^)(oo). However, in the KA models it is possible that 
some particles will never be able to move but are still 
culled because other particles that may move but block 
them are culled, and thus Al¬ 
though the case P^^^^\oo) > = 0 is possible 

in finite systems, we assume that it does not occur in the 
thermodynamic limit since we encountered such a seen 
ario only in very small systems. 
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III. CULLING DYNAMICS 

A. Culling Dynamics for the m — 1 FA and KA 
models 


In the m = 1 models, there are no permanently frozen 
particles, and thus M (oo) = 0. Furthermore, we obtain 
an explicit expression for M{s). The number of particles 
culled in the s’th step, M{s — 1) — M{s), is the number 
of particles that all of their {s — l)-nearest neighbors are 
occupied and at least one of the s-nearest neighbors is 
vacant. For d = 1 this is 


- 1) - (1 - p^) , (4) 

and for d = 2 it is 

(s - 1) - (s) = pMs-i) (1 - /"). (5) 

Solving these recursion equations yields for d = 1 

(6) 

and for d = 2 

M(24)(s) = p25(s+l)^ (-7^ 

Hence, by substitution in Eq. (I3| we find that the culling 
times are given by 


r(bi) = 


l-p2’ 


( 2 , 1 ) ^ ^p2s(s+l) 


s=0 


02 (0,P^) 


(8) 


where 02 is the Jacobi Theta function [55|. At high 

(2 1 ) 

particle densities, <C 1, we may approximate Tc ’ by 


ri-' . {I- (9) 

In a similar manner for general dimensions, by cal¬ 
culating the number of 5-nearest neighbors in a d- 
dimensional hypercubic lattice, Gd{s)^ we find that 

^(d.l)(g) (10) 

It was shown in [s^l that Gci(s) is a polynomial given by 

= (a) 

k=o ^ \ / 

At high particle densities, 'r’ = l — p^l, we can find 
an approximation for Tc in any dimension. From Eq. m 
we find that is given by 


s=0 s=0 


We now note that Gd{s) is a polynomial of order d with 
the coefficient of given by 

G,is) ='t{t)r ‘‘+^ 

k=0 ^ • 

(13) 

Changing the sum over s to an integral over x = 
yields 


poc 

^(d,l) ^ y-l/d / 

Jx=0 

_ r (1 +1) (dpv-^ 

2i;Vd 


(2^ 

d\ 


dx = 


(14) 


Except for the non-trivial prefactor, the dependence 
of Tc on the vacancy density v comes simply from the 
fact that Tc is the distance to the nearest vacancy, which 
scales as 


B. Culling Dynamics for the m = 2 FA and KA 
models 

In the d = m = 2 models, we find and 

(2 2 ) 

Tc ’ numerically by running simulations on square sys¬ 
tems of size L X with L = 100 or 1000. We only 
consider densities below the critical density (see Eq. 
(dD), pc{L = 100) ^ 0.94 and pc{L = 1000) 0.96, 

since in the thermodynamic limit the critical density is 
Pc {L = oc) = 1 and thus the results relevant to this limit 
are below the size-dependent critical density. 

Eigured] shows the dependence of M(s) on s. At small 
s we find an exponential decay — InM^^’^^ ~ <s, which is 
similar to the behavior of in Eq. (j6j), while for 

large 5, the form is Gaussian — InM^^’^^ ^ 5^, which is 
similar to in Eq. 0. The reason is that for small 

s, the particles are culled mostly one by one such that the 
behavior is quasi-one-dimensional, and when the empty 
region is large enough, the particles around it are culled 
by diagonal shells as a two-dimensional system, see Eig. 

m 

In the m = 1 models, the culling process is equiva¬ 
lent to expanding the vacant regions, and thus the culling 
time for a given site is its distance to the nearest vacancy, 
Tc see Eq. m above. In the m = 2 models, 

the culling process is dominated by critical droplets [52| , 
which are small regions that may be expanded by the 
culling process to include the entire system. Hence, the 
culling time is the distance to the nearest seed for a crit¬ 
ical droplet, if the site is far enough from the droplet, as 
shown in Eig. [3l Because the probability of a given site 
to seed a critical droplet is exp (—2 A/t), with A ^ 0.25 
in the sampled density range Q, the average distance 
from a droplet, and thus the mean culling time, should 
scale as exp (A/t). We see from Eig. |4]that this form of 
the scaling is consistent with our numerical data. 
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FIG. 1: The fraction of sites that haven’t yet been culled, M, 
(top panels) and its logarithm (bottom panels) as a function of 
the iteration number s for the d = m = 2 FA model (left pan¬ 
els) and KA model (right panels). The system’s linear size is 
either L = 100 (blue dots) or L = 1000 (red continuous lines). 
Each curve is for a different density (from left to right) p — 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.91, 0.92, 0.93, 0.94, 0.95. 
There is almost no difference between the different system 
sizes. The bottom panels show only data for L — 1000. The 
straight black lines in the bottom panels are ~ s and ~ . 

There is almost no difference between the FA and the KA 
models. 


FIG. 3: A scatter plot of the culling time s vs. the dis¬ 
tance (I Ax I + |Ay|) to the nearest critical droplet for a sin¬ 
gle 100 X 100 configuration at the critical density p — 0.951 
for the d=m=2 FA model. Lines with positive (negative) 
slope indicate culling around the empty region away (toward) 
the seed of the droplets, see text. At intermediate times 
(50 < s < 150) the culling time generally grows with the 
distance from the droplet, since at those times the droplet in 
the configuration shown here expanded away from the seed 
in two of the four directions. At longer times, the droplet 
expanded in all directions, including toward the seed. The 
points with low s and large distance represent small, active 
regions which are far from the critical droplet. 



FIG. 2: The culling time s for each site in a typical 100 x 100 
configuration in the FA model at p = 0.92. The legend shows 
the range of s represented by each color. The diagonal borders 
between regions of different colors, which are indicative of 
a two-dimensional culling process, are clearly seen at large 
scales. See for example the diagonals at the lower left corner, 
and in the region above the legend to the right. 



FIG. 4: The mean culling time Tc vs. the reciprocal of the 
vacancy density l/xford = m = 2. There is almost no dif¬ 
ference between the KA and the FA models. Below the jam¬ 
ming density, there is almost no difference between L = 100 
and L — 1000, and the mean culling time may be fitted to 
Tc = Aexp(A/x) with A ^ 2.1(KA),1.5(FA), and A = 0.25, 
as shown by the black dotted lines. Tc reaches a maximum 
around the critical density (vertical dashed black lines), be¬ 
cause at higher densities the frozen sites can no longer be 
culled and thus do not contribute to Tc. 
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IV. PERSISTENCE IN THE PHYSICAL 
DYNAMICS 


A. Real Dynamics for the m = 1 FA and KA models 


Since the m = 1 KA lattice gas is equivalent to SSEP, 
instead of considering motion of particles, we may think 
of the dynamics as diffusion of vacancies, such that the 
persistence of a given site is the mean first passage time 
of vacancies to that site. Since the vacancies can dif¬ 
fuse freely, all particles will eventually move and thus 
P(oo) = 0. At high particle density, ^ 1, we may 
make the approximation that the vacancies are indepen¬ 
dent, and allow two (or more) vacancies to occupy the 
same site. Furthermore, at long times the discrete na¬ 
ture of the lattice becomes irrelevant and we may use 
results from continuous models. Under these approxima¬ 
tions, the long time behavior of the mean first passage 
time distribution, and thus of the persistence function, is 
described by [Hsl-I^ 




exp 


exp 


-2vx ^ 


AirDvt 
■ In Dt/R^ 


d = l 

d=2 > (15) 


exp [— {d — 2) SdR^ ‘^Dvt] d > 3 


where Sd = jV {dl2) is the surface area of the d- 

dimensional unit sphere, and D is the self-diffusion coef¬ 
ficient for the motion of the vacancies. For the m = 1 
models, D = l/(2d). In continuous models, R is the ra¬ 
dius of the trapping region. In a discrete lattice, in which 
each site contains at most one particle, P = 0(1). 

At short times we may use a mean-field approxima¬ 
tion, such that the probability that a particle can move 
to an adjacent site (for the first time, since this is an 
approximation for short times) is and thus 

dP 

= -vP{t), (16) 

which yields P{t) = Note that in kinetically con¬ 

strained models, by construction the occupation proba¬ 
bilities of neighboring sites at a given time are uncorre¬ 
lated, the dynamics are spatially heterogeneous [tI. [3^. 
In the appendix we derive an exact expression for the 
one-dimensional case at all times, under the approxima¬ 
tion that the diffusing vacancies are independent. 

The FA Ising model may be thought of as a diffusion- 
reaction model, which behaves similarly. At high parti¬ 
cle densities in the FA model, lUio ^ lUoi, and sites in 
state 0 can be considered to change practically instantly 
to state 1 (when the kinetic constraint does not prevent 
them from doing so) compared to the time it takes a site 
in state 1 to flip. Thus, when a state 1 flips and imme¬ 
diately after that its state 0 neighbor flips, it appears as 
if the state 1 moved. This effective movement happens 
on a different time scale than in the KA model, because 


the rate lUio is smaller than 1, and thus time should be 
normalized by WoiIUio = min m order to map 

the FA dynamics on those of the KA model. In what 
follows we thus normalize time by lUoiIUio and interpret 
the rates lUio and lUoi as equal to unity in the KA model. 

Figure [5] shows the persistence function for the m = 1 
models. We see that the FA and KA models behave 
similarly, except for a prefactor, and that the analytical 
approximation, Eq. m, is in good agreement with the 
numerical results. 
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FIG. 5: The persistence function P as a function of the nor¬ 
malized time WqiWiqI in the m = 1 models in d = 1 (a) and 
d = 2 (b) for the FA (dashed lines) and KA (dotted lines) 
models. In panel (a), the continuous lines are the analytical 
approximation, Eq. ([T5l) . As the scaling from Eq. ([T5l) for 
d = 2 is difficult to see from panel (b), we show in panel (c) 
t/ InP vs. t ior p — 0.9999 in the d = 2 FA model. At small 
times P decreases exponentially with t, but at longer times 
it decreases exponentially with t/lnt as expected from Eq. 
(|15l) . The drop at long times is due to finite-size effects, and 
the thin line is ~ Int. 
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B. Real Dynamics for the m — 2 FA and KA models 

Similarly to the m = 1 models, at high densities the 
m = 2 FA model behaves as the m = 2 KA model under 
the proper time normalization for the exact same reasons 
as in the m = 1 models. At short times, the persistence 
decays exponentially as 

P{t) = (17) 

where in d = 2, 

(18) 

is the probability that a random particle has at least two 
neighboring vacancies before and after the move. This 
result is obtained from a mean-field approximation, valid 
for short times, similarly to the analysis above for the 
m = 1 models. This exponential decay continues until 
time 1/ITio, which is the typical time at which all the 
sites that were able to flip at t = 0 have flipped (FA) or all 
the particles that were able to move at t = 0 have moved 
(KA). After that time, we see from the numerical results 
that the persistence decays as a stretched exponential 


_Ln(p(2.2)) 




P(t) = 


(19) 


as shown in Fig. [6l Note however that in the thermody¬ 
namic limit and at extremely long times the persistence 
function eventually decays exponentially ji^ . 

Combining Eqs. (0), dni) and m we find that the 
persistence time may be approximated by 


WioTp = 


l-e-“ 

n 


r(i ii) 


( 20 ) 


where F (a, z) is the incomplete Gamma function. In the 
limit of small v (and thus small /3), Tp may be further 
approximated by 


I O / -I \ 1//5 

Previous numerical studies 0 have shown that in a 
50 X 50 system, the exponent [3 converges to a value of 
0.42 as the density is raised to p = 0.917, near the critical 
density for jamming in a system of that size, pc ~ 0.94 
[63|, while we get /3 = 0.25 at p = 0.92 and it clearly does 
not converge. The reason for this apparent discrepancy 
lies in the preparation protocol. In our simulations, each 
site at time t = 0 is in state 1 with probability p and in 
state 0 with probability and thus the initial configura¬ 
tion is chosen from the equilibrium distribution. In the 
simulations reported in [^, all the spins were initially 
set to 0, the system was evolved for a long time until it 
apparently reached equilibrium, and then the measure¬ 
ment of the persistence function started. However, we 
suspect that these simulations did not equilibrate. In¬ 
deed, by simulating a 50 x 50 system at p = 0.92 with 


FIG. 6: (a) The persistence function P as a function of the 
normalized time in the m — 2 FA (dashed lines) and KA 
(dotted lines) models. At short times the exponential decay of 
the persistence function agrees with Eqs. GniHj, denoted by 
the solid lines, but at longer times it decays sub exponentially 
as —InP ~ (b) The exponent (3 vs. the vacancy density v. 

The blue symbols are the numerical results, and the dashed 
red line is the analytical approximation for small Eq. (1361) . 


such a quenched initial condition, we find yd = 0.42 if we 
wait for 10^ steps per site, while by waiting for 10^ steps 
per site we get /3 = 0.32, still far from the equilibrium 
result. It would be interesting to test the convergence of 
yd to its equilibrium value by waiting considerably longer 
times after such quenches. 


V. COMPARISON BETWEEN THE CULLING 
DYNAMICS AND THE REAL DYNAMICS 

In order to compare the culling time Tc and the per¬ 
sistence time Tp, we return to the picture of diffusing 
vacancies. Eor m = 1, combining Eqs. and m 

for d = 1 and d > 3 yields 


Ml) _ 

Ip 


TTT^ d = 1 
CdTc d>3 

where q is some constant. Eor d = 2 the integral 


( 22 ) 


,-( 2 , 1 ) = 

Ip 


/ oo 

p(2-l)(i)dt, 


(23) 
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where is given by Eq. m, cannot be computed 

exactly for any finite t’, and requires more work to find 
the asymptotic expansion for small v. We first change 
the integration variable from t to x = AirDvt 


r(2,l) 


-/ 

AttDv 74, 


exp 


AttDv 


\l].{x/c2v) 


dx, (24) 


Solving for v yields 


4/3A 


{P - 3) W 


/ 3-3 


(576eA3/3(5£,) 


/3/2\ V(3-^)' 


(33) 


where 


C2 = 47rP^, 


where W [^] is the product-log function [65| defined as 
(25) the solution to 


and as noted above P = 0(1). We now divide the range 
of integration to three parts: AttDv to 1, 1 to l/(c 2 'r’), 
and l/(c 2 'r’) to 00 . The first part is negligible because 
its total contribution is smaller than 1. The third part 
is negligible because in the limit of v 0 , both the 
integrand and the range of integration go to 0. Hence, 


^(2,1) 

Ip 


I r^lic2v) 


X 

Inx — lii{c2v) 


dx. (26) 


In this region, |lnx| < |lnc 2 r’|, and thus we may further 
approximate by 


.( 2 , 1 ) 


AttDv 

ln(c2i;) 

AttDv 


1/{C2V) 


exp 


ln(c2i;)_ 


dx 


(27) 


where in the last approximation we used v 1. Com¬ 
bining Eq. ([TT|) and ([27|) yields 


ryi)~r 2 lnre. (28) 

In the m = 2 models, the dynamics is dominated by 
movement of droplets, not of individual vacancies. We 
recall that the droplets appear with an effective density 
of 


h = exp (—2A/'z;), (29) 


z = W [z] exp {W [z]) . 

Expanding Eq. (ESI) for small (3 yields 

3A/3 

21n(l//3)' 


(34) 


(35) 


Solving for /3, and approximating for small v yields 


/3=^ln 

^ 3A 



(36) 


which for small v is shown on Eig. [ 6)3 to roughly agree 
with the value of /3 obtained from the stretched exponen¬ 
tial form of the persistence function, Eq. m- 

In order to show that the relation between Tp and Tc is 
general, we also consider here the two-dimensional spiral 
model [l3, [l5|. This model jams at a finite density, at 
which the frozen structures are one-dimensional strings 
that run along the diagonal directions of the lattice. Be¬ 
low the critical density, but near it, the largest contri¬ 
bution to the persistence time comes from such almost- 
frozen strings, and thus we may approximate this as a 
quasi-one-dimensional process (see Eig. Q, which leads 
to 


rsvival ^^2 jd, 


(37) 


and that ^ exp (A/i;) ^ v The self-diffusion coef¬ 
ficient of particles in the m = 2 model is given by Q 

D = exp {—2X/v) ^ (30) 

We are interested in the self diffusion of droplets. The 
particles inside the droplets are the most mobile parti¬ 
cles in the system, and thus contribute the most to the 
self diffusion coefficient of particles in the system. There¬ 
fore, the self diffusion coefficient for the droplets may be 
approximated by the self diffusion coefficient of the par¬ 
ticles, given by Eq. EOl)- Using Eq. (|3Q]) and changing v 
to V in Eq. dzni yields 

4^’^^ "Tc lUTc. (31) 

Erom this relation we can find an approximation for 
/3{v) at small v. Combining Eqs. ([ 21 ]), (j30|) and (|3T]) for 
the KA model yields 


As the density approaches the critical density, we see 
from the numerical results shown in Eig. El that the 
diffusion coefficient approaches zero as 

Dr^r-\ (38) 

and thus 

^spiral _ ^ 4 _ ( 3 g) 

Figure El shows the excellent agreement between the 
numerical results and Eqs. El), (EHl, (EB) and (|39|) at 
densities slightly below the critical density {pc = 1 for the 
EA and KA models, and pc ~ 0.7 for the spiral model). It 
would be interesting to study the structural connections 
between the directed percolation underlying the jammed 
structures in the spiral model and the one-dimensional 
nature of the relaxation processes in it. Eurthermore, 
it would be interesting to identify possible logarithmic 
corrections to Eq. (|39|) and to numerically test their ap¬ 
plicability. 
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FIG. 8: The self-diffusion coefficient D in the two-dimensional 
spiral model vs. the cuffing time Tc. The continuous line is 



FIG. 7: Spatial structure of the cuffing (top) and persistence 
(bottom) times in the spiral model at p = 0.61, which is 
slightly below the critical density for the simulated system 
size of L — 100. The one-dimensional structures running 
along the diagonal directions of the lattice are clearly seen. 


VI. SUMMARY 

In this paper we investigated the relation between the 
structural changes in the system, represented by the 
culling time Tc, and between the relaxation time of the 
persistence function in the Fredrickson-Andersen and 
Kob-Andersen kinetically constrained models in one and 
two dimensions, and in the two-dimensional spiral model. 
We found that ^ rj, up to logarithmic corrections in 
the Fredrickson-Andersen and Kob-Andersen models in 
two dimensions, where 7 is model-dependent. This re¬ 
sult is explained by mapping the persistence of a site to 
a first passage time of diffusing defects, with their initial 
distance given by the culling time. 

We also found that the persistence function in the 


WoiWioTp 



FIG. 9: The relation between the normalized persistence time 
WioWoiTp and the culling time Tc for the FA and KA mod¬ 
els and for the spiral model. The continuous lines are the 
analytical approximations, Eqs. (Eai, (HU, (EU, and (j39p . 


m = 2 models at long times behaves as a stretched ex¬ 


ponential exp 


{t/tof 


where [3 probably goes to zero 


at small vacancy densities in contradiction to previous 
studies in which [3 was believed to converge to a finite 
value. The difference arises because the previous results 
were obtained in systems which were not in equilibrium, 
while our simulations are performed in equilibrium. 

The general relation between the culling and the per¬ 
sistence times may also hold in other models, including 
continuum models, and in experiments. Since the m = 1 
Fredrickson-Andersen and Kob-Andersen models repre¬ 
sent normal gas or liquid, while the m = 2 models and 
the spiral model represent glassy behavior, the exponent 
7 is a measure for the “glassiness” of a system. It would 
be interesting to check whether this relation holds also for 
other measures of the relaxation, such as the autocorre¬ 
lation function, with the same value of 7 . Also, it would 
be interesting to study the relation between the culling 
and the persistence in the three-dimensional extension 



















of the spiral model in which there is a decoupling 
between the structure and the dynamics, namely the den¬ 
sity above which permanently jammed structures appear 
is lower than the density above which the long time self¬ 
diffusion stops. 
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where In{t) is the modified Bessel function of the first 
kind. Setting the general solution in the equation for 
k = 1 yields 

oo 

^ Aih-kit) = 0 . (A5) 

Z= —OO 

Using the relation In{t) = we find that 
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Appendix A: Exact expression for P{t) in the 
one-dimensional m — 1 KA model 


Here we derive an exact expression for the persistence 
function in the one-dimensional m = 1 KA model, using 
the approximation of non-interacting diffusing vacancies. 
Consider a one-dimensional lattice of length 2L -\-l with 
V • (21/ + 1) vacancies. We are interested in the limit 
1/ ^ oo, and implicitly take this limit during the deriva¬ 
tion whenever there is no singularity. The persistence 
function, P(t), is the probability that none of these dif¬ 
fusing vacancies reached the origin until time t. Since 
we assume that the vacancies are non-interacting, it is 
enough to compute the probability that a single vacancy 
did not reach the origin until time t, Q(t), and from that 
we can obtain P(t) = 

At each time step S, there is a probability S that the va¬ 
cancy tries to move, and if it does there is an equal prob¬ 
ability to move either to the left or to the right. Without 
loss of generality, we may assume that this vacancy is at 
site /co > 0 at time t = 0. The evolution equation for the 
probability Vk (t) that at time t the vacancy was at site k 
reads 

Vk{t + S) = (1 - 6)vk{t) + - [vk+i{t) + Vk-i{t )], (Al) 
for k > 2, and 

Vi{t + S) = {1 - S)vi{t) +-V 2 {t), (A2) 

for /c = 1, since if the vacancy reached the site /c = 0, the 
process stops. In the limit ^ > 0 this transforms to the 

differential equations 


dvk 1 . , . 

= -Vk + - + Vk-l) , 

dvi 1 


(A3) 


The general solution to the first differential equation is 


oo 

^ (A4) 

l = — C)0 


Al = (A 6 ) 

Imposing the initial condition r’/c(0) = 6k,ko and using 
^n(O) = ^n,o yields 

Ak = 6k,ko - 6k,-ko‘ (A7) 


Hence, 


Vk{t) = e * [Ik-ko{t) - Ik+ko{t)] ■ (A 8 ) 


The probability that the vacancy did not reach the 
origin until time t, given that it started from ko is 


Qi{ka,t) = '^Vk{t) = I 


*=1 


ko — 1 

Iko (t) + ^0 (t) + 2 Ik (t) 
k=l 

(A9) 


Averaging over all initial states yields 



1 X ^ 

- Qi{ko,t) = 

ko = l 

-j _^ Y f^ kQ 1 

2- ^- ^ ^/o(t) + 2e ^ Yj , 

ko = l k=l 

(AlO) 


where we used 


(All) 

n=l 

In order to calculate the last sum, we change the order 
of summation such that 


L ko — l oo L 

E E 'M = E E '‘(*) = 

ko = l k=l k=l ko=k-\-l 

oo y- oo 

= Y,{L-k) hit) = 2 (e* - k) - (A12) 

k=l k=l 

where we used Eq. (jAlip . 

We now use the relation 

kikit) = ^ [Ik-i{t) - Ik+i{t)], (A13) 

and find that 

oo 

Y,kIk{t) = -[Io{t) + h{t)]- (A14) 

k=l 
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Therefore, by combining Eqs. (ITTnli . (ITT^ and (IA14II . 
we find that 


Q{t) = ^ ^ ® - te * {Io{t) + /i(t))] , 

(A15) 


and thus the persistence function is given by 
P{t) = lim = 

L^oo 

= exp{-v [e“*Jo(t) - 1 + te~* {Io{t) + /i(t))]} . 

(A16) 


For t ^ 1 the persistence function behaves as 

exp (—ut), and for t ^ 1 it behaves as 


P{t > 1) 



(A17) 
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